Why is d — wave pairing in HTS robust in the presence of impurities? 
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In the recent theory of strong correlations by |l| - it was shown that by lowering doping 
concentration a forward peak is developed in the charge scattering channel. Accordingly, near the 
optimal doping the nonmagnetic scattering is pronounced in the d-channel and d — wave pairing 
is robust against defects and impurities. It is reflected in a decrease of the slope of the critical 
temperature Tc{Ts) at Fs ^ by factor (1 — /3), where f3 is the scattering anisotropy parameter and 
Ts is an average scattering rate. For large doping, /3 is small Q , and d — wave pairing loses its 
robustness. The theory is generally formulated for the bi-layer model by including: (1) intra- and 
inter-plane pairing; (2) intra- and inter-plane impurities. 



Recently there were a lot of evidences for d — wave 
pairing in high-Tc-superconductors (HTS). However, the 
mechanism of pairing is still unknown, although many 
proposals were done. For instance, there were claims 
that superconductivity is due to strong correlations [Q, 
or that AF fluctuations are important - see Even 
the electron-phonon interaction renormalized by strong 
correlations can give rise to d — wave pairing, as it was 
shown by Kulic and Zeyher Moreover, in 

0, 1^, 1^ it was shown that the charge-charge scatter- 
ing is renormalized via a charge- vertex 7 ((7) (screening 
due to strong correlations) with a peak in the forward 
scattering (| q |^ kp), while the backward scattering 
(large | q \) is suppressed. This means that a bare non- 
magnetic scattering potential Uo{q) is renormalized by 
strong correlations Uo{q) Ueff{q)- Note, in the Born 
approximation (after averaging over impurity positions) 
one has Ueff{q) = ^^{q)u'^{q) = u'^iq), where uo{q) is 
the single impurity potential. It follows 0], 1^ that 
when the bare potential is local, i.e. Mg(g) — const, Ueff 
becomes " long-ranged" in the presence of strong correla- 
tions. Moreover, within the t~J model |^ and at 
some optimal doping (5 « 0.2 the d-channel component 
of Uef f (q) is of the same order of magnitude as the one 
in the fully symmetric s-channel. When Uo{q) is due to 
the electron-phonon interaction then superconductivity 
can be also of d-type in the presence of strong Hubbard 
repulsion. 

If we accept that superconductivity in HTS is of 
d — wave type, then one can put a question: why 
is d — wave pairing robust in the presence of various 
kinds impurities and defects in HTS materials, as it 
was observed experimentally (|], 0? In this paper we 
present a simple theory, which is based on the effec- 
tive nonmagnetic impurity potential Ueff{q) with the 
large d-c_omponent and investigate its effect on d — wave 
pairing M. Experiments on the irradiation of samples as 



well as on substitution of atoms in HTS materials are 
analyzed by using the proposed theory. In the following 
we consider the weak coupling theory |9| and its exten- 
sion to the bi-layer systems [0, [0. In the bi-layer 
model (adequate for YBa2Cu'iO'j-x) with the hopping 
t±_ between two layers (separated by Y ions) and with 
the intra- and inter-layer pairings A||(p) and Aj^(p) re- 
spectively the averaged (over impurity positions) Green's 
function G{p){= G{p,iujn)) is 4 x 4 matrix. It is given 
by G^^ip) = Gq ^(p) - %mp{p)G{p), where G'o(p) is the 
bare Green's function and the impurity self-energy ma- 
trix Y,imp{p) in the bi-layer model is given by 
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lJ:jf{p,^)G{pf,lOJn). (1) 



Here, a;„ — 7rT(2n -I- 1) and c are Matsubara frequencies 
and the concentration of impurities respectively. In or- 
der to simphfy notations we set UlJ'^{p,p') = U{p,p'). 

The structure of the 4x4 matrix U{p, ff) depends on the 
type of impurities: Case / - there are inter-layer impu- 
rities (for instance substitution of Y by Pr); Case // - 
there are intra-layer impurities (for instance substitution 
of Cm by Zn, or the presence of oxygen vacancies due to 
ion damaging etc.). They are explicitly given in [ p^ - see 
also below. Since, in the weak coupling limit the rele- 
vant processes are taking place near the Fermi surface it 
is assumed that the matrix U{p,p') = U{d,d') depends 
on the angles 9, 0' of momenta on the Fermi surface. In 
the following, the definition J^^ = {(■■■))e' is used. 

G~^{p, iujn) in the bi-layer model is searched in the form 
analogous to that given in with renormalized 

Matsubara frequencies ujn and renormalized supercon- 
ducting order parameters A„^i (p) {i — 1,2 is the band 
index -see below) 



W„(0) =CJ„ + i(Tj/( 



9') 



-)e' (2) 



(3) 

To this set of equations one should add also the self- 
consistent gap equations 

AM = TcJ2{K](.0,e')^^M^^=)e'. (4) 

V ^uUO') + Alj{e') 

The summation over the bar indices is assumed every- 
where. Here, Aj(6') {i = 1,2) is the gap in the i — th 
band, i.e. Ai^zie*) = A|| (p) T Aj.(p) where X^J{9,0') 
(i, j — 1,2) are the intra- and inter-band pairing coupling 
constants |p^ . In the following, we assume d~wave pair- 
ing (instead of the s — wave one used in |l^) Ai{0) = 
A,Yd{e) = -A,(0 + V2), {Yd{e))e - 0, {Yi{0))e - 1, 
and \,j{e,e') = KjYd{6)Yd{0'). Near the critical tem- 
perature Tc one can linearize eqs.(2)-(4) which gives 

(i„ (0) = C^„ [1 + - ^ " / ' ] = UnVn {e) (5) 



n ' ' 



where ?7„ = l-l-(rs4i-f rs^22)/ | |, r?n,ij = ^»i+rdjj/(| 

I -t-r^-Td), = rs,ii-i-rs^22, and = rd,ii-i-rd,22- 

Eq.{7) contains reach physics which depends on the val- 
ues and signs of the coupling constants Ay see . How- 
ever, in order to illustrate the basic effect of strong cor- 
relations we simplify the problem by assuming that the 
inter-layer pairing is negligible, i.e. A_i(p) « and 
Ai « A2 ~ A||. Then eg. (7) is reduced to 



in^ = ml)~m-P>+l)], (8) 



^nM = Vn,,j{0)Aj{0), (6) 

where T~j^{0,0') and ?7„,ij(^) depend on the type of im- 
purities (intra- or inter-plane) -see |l^, while ri„ ij{0) are 
solutions of a coupled set of integral equations (w.r.t. 
0). In the presence of inter-plane impurities (the case 
/), one has t-/{0,0') ^~J~N~(0ju{0,0'), where N^{0) 
are the fully symmetric angle dependent band (i = 1,2) 
density of states on the Fermi surface. Note, this kind 
of impurities does not affect isotropic s — wave pair- 
ing - see ]l^, but it does affect d — wave pairing [p^ . 
In the case // of the intra-plane impurities one has 
T-/{0,0') = Ui{0,0')Nj{0) with Ui{0,0') = U2{0,0'). 
Note, the s — wave inter-layer pairing is strongly affected 
by these impurities, due to the breaking of the local re- 
flection symmetry between layers - transitions be- 
tween symmetric and antisymmetric band occur. 

Let us analyze the case // . which is frequently realized 
in the YBa2Cu307-x family. A very important ingredi- 
ent of our theory are strong correlations, for which it has 
been shown that in t~^(0, 0') only s- and d-channels dom- 
inate 1^] , , Q . When generalized to the bi-layer model 
one has Tr^\0,0') - T,^,,[0,0') + Td^,,Yd{0)Yd{0') + .... 
Note, T s^ij{0T0'){> 0) is fully symmetric and it holds 
Ts,^,{0,0'') > Td,^JYdi0)Yd{0'), where Td^y > 0. Their 
doping dependence is analyzed below. Let us assume 
that the fully symmetric scattering rates are weakly an- 
gle dependent, i.e. Ts^ij{0,0') ~ Ts,ij- In this case eq.(3) 
is reduced to a set of algebraic equations for the critical 
temperature Tc 



where (3 = (Fd/F^) < 1, x = TjAnTc and T° = 
l.lSwcexpj— 1/A}, A = All + '^12- Similar equation 
has been obtained also in |l3j. It tells us that the 
pair-breaking effect of nonmagnetic impurities is strongly 
reduced if the scattering potential contains a large 
anisotropic part which is compatible with the pairing am- 
plitude. It follows from eq.(8) that in the presence of an 
anisotropic scattering {P ^ ) d — wave pairing is more 
robust than in the case of the isotropic one (/? = ). 
This robustness is reflected in the decrease of the slope 
of Tc(x) at X = by factor (1 — f3). It should be stressed 
that the order parameter keeps its d — wave form, i.e. 
A{0) « AdYdi0) and {A{0)) = 0. 

The above simple analysis would not deserve to much 
attention if we would not be able to find the range of 
the parameter as well as the origin for its large value 
in HTS. Luckily, the doping dependence of /3(5) can be 
extracted from the theory of strong correlations jl], 
Q, where for 6 < 0.2 the parameter /? becomes of the 
order of one(!), while at large doping {1 — S 0) one 
has /3 <^ 1 -see Fig.l. The latter value leads to a rapid 
decrease of Tc by increasing F^. These results tell us that 
d — wave pairing in the underdoped and optimally doped 
HTS materials is robust in the presence of impurities, 
while it is much less robust in the overdoped systems - 
see Fig.l. 
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see Fig. 2. 




FIG. 1. The anisotropic scattering parameter (3 as a func- 
tion of the doping concentration S based on the work . 

The average pair-breaking parameter is extracted 
from the measurements of the residual resistivity pi — 
An/TtrUJpi, which is due to the substitutional defects 
{Cu — > Zn, Y — >Pr) or due to the ion (A^e+) or electron 
damaging of HTS materials - see [Q. The plasmon 
frequency cupi is extracted from the slope a in the linear 
resistivity p{T) = pi + aT, where {An/uj^i) — a/2n\tr 
and Xtr is the transport coupling constant. The realistic 
guess for Xtr is Xtr ~ 0.3 However, it is not correct 
to replace Tg by l/rtr because for the latter backward 
processes dominate. It turns out that in the presence 
of strong correlations one has Fg = p • l/rtr, i-e. Fs 
= 2pT:XtriPi/a) and for/below optimal doping S < 0.2 it 
is obtained p « 2 and /3 > 0.8 |], §. 

We discuss here three scenarios for the suppression of 
the measured slope {dTc/dx)x=a {x — Pi/a), which are 
compared to the prediction of the A — G theory |^] 
for d — wave pairing with isotropic impurity scattering. 

1) d — wave pairing {{A{9))g = 0) and anisotropic 
scattering {(3 ^ 0) 

Based on the above analysis and from eq.(8) one has 



{dTjdx)x=o = - /3) V^VS, 



(9) 



where x — Pi/a. The experimental value is 
I dTjdx U^o« 0.25 @, §, 0, |18). For p « 2 
and Xtr ~ 0.3 ^ one obtains « 0.85. Such a large 
experimental value of (3 is in accordance with the predic- 
tion of the theory of strongly correlated electrons 0,0, 
. This means that the present theory can consistently 
explain the experimental slope | dT^/dx U=o 111, - 
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FIG. 2. The critical temperature Tc[K] of d — wave su- 
perconductor as a function of the experimental parameter 
Pi/a[K]. The case /3 = is the prediction of the standard 
d — wave theory with isotropic scattering. The experi- 
mental data are given by crosses - YBa2Cu3_a;Zny07_5 
and circles - Yi.j^Prj^BaaCuaOT-a H]. 

2) Anisotropic pairing ((A(6'))e ^ 0) and 
isotropic scattering {(3 — 0) 

Based on eqs.(2)-(5) one obtains eq.(9), where (3 — > 
[3an = (A(0))2/(A2(0)) - scc 0. For I dTj dx U^o^ 
0.25 one has (311 « 0.85, i.e. (A(6'))e « {l\'^{e)). 
This result means that a large s — wave component ap- 
pears in the order parameter, i.e. A(0) = A^ -I- /S.dYd{d) 
and Ad ~ 0.5As. This result is compatible with the 
experimental observation of the finite critical current in 
the Ph — YBCO Josephson contact along the c — axis 
]l6t , but it is incompatible with d — wave like pairing 
(I {A{0))g I /Aniax ^ 1) fouud in many experiments. 

3) Anisotropic pairing {{A{6))g ^ 0) and 
anisotropic scattering (/? ^ 0) 

In this case one should set /3 (3eff in eq.(9), where 
/3e// = ml + AD /jAl + Al ). Since Plf^ « 0.85 one 

obtains A^ = 2Ad\/l — (3/0.85. This means that finite f3 
suppresses the s — wave component in A{9). In our opin- 
ion, the realistic situation for the optimal doping might 
be described by the case 1), or by 3) with 0.78 < (3 < 0.85 
and < As < Ad/3. 

In conclusion, the effect of nonmagnetic impurities on 
d — wave superconductivity in HTS materials is stud- 
ied within the theory of strong correlations [|l| , ||^ , |^ 
which predicts: (i) by lowering doping d the forward 
scattering peak is developed in a nonmagnetic scatter- 
ing where a large d — wave component is present in it; 
(a) transport properties are renormalized in such a way 
that l/rtr = l^s/p, where p > 2 for 5 < 0.2. Based on this 
theory, we have shown that the robustness of Tc in HTS 
systems in the presence of defects and impurities is due 
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to: (a) an anisotropic scattering (due to strong correla- 
tions) and (b) the anisotropic pairing with {A{9))g = 
or I {A{0))e I /Amax 1. The theory, which is gen- 
erally formulated for the bi-layer model with intra- and 
inter-laycr pairing and with intra- or intcr-laycr impu- 
rities, shows that the underdoped and optimally doped 
systems are robust in the presence of impurities, while 
the overdoped systems lose this robustness. 
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